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The motion of light and a neutral test particle around the charged D-star has been studied. The 
difference of the deficit angle of light from the case in asymptotically flat spacetime is in a factor 
(1 — e 2 ). The motion of a test particle is affected by the deficit angle and the charge. Through 
the phase analysis, we prove the existence of the periodic solution to the equation of motion and 
the effect of the deficit angle and the charge to the critical point and its type. We also give the 
conditions under which the critical point is a stable center and an unstable saddle point. 
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(~~) . 1. Introduction 

(N ; 

w , Astronomical observations on the CMB anisotropy[l-3] and the relation between red-shift and luminosity distance 
of SNela [4-6] depicted that our universe is spatially flat, and about 23 percent and 73 percent of energy density is 
U I resulted respectively from cold dark matter and dark energy, respectively. The nature of these substances are quite 

■ unusual and there is no justification for assuming that it resembles known forms of matter or energy. On the other 
hand, phase transitions of quantum fields in the early Universe may produce various kinds of topological defects. The 
idea that point-like defects known as monopoles ought to exist has proved to be remarkably durable. The global 
monopole, which has divergent mass in flat space is one of the most important defects. The effects of gravity on the 
global monopole were firstly studied by Barriola and Vilenkin [7]. When gravity is taken into account the linearly 
divergent mass of a global monopole has an effect analogous to that of a deficit solid angle plus that of a tiny mass 
at the origin. It is shown that this small gravitational potential is actually repulsive [8-12]. 

As possible candidates for the cold dark matter, various kinds of cold stars such as Q-stars have been proposed [13- 
22]. Later, a new class of cold stars named D-stars (effect stars) have been proposed by Li et al. [23,24]. Compared to 
Q-stars , the D-stars have a feature, that is, in the absence of the gravitational field the theory has monopole solutions 
""q ' , which makes the D-stars behave very differently from the Q-stars. Recently , these new objects have farther been 
qh, studied. Their possible interesting astrophysical applications [25,26] and charged D-stars have also been studied. In 
^ 1 contrast to ordinary D-stars, there exists a U(l) gauge field outside the charged D-stars [27]. 

In this paper, we will discuss the behavior of test particle and geodesies in the exterior spacetime of the charged 
I* ■ D-stars. With the aid of phase-plane analysis, we carefully analyse the motion of light and a test particle around a 
. ' charged D-stars and show some unique features of the behavior. 

x : 

2. The charged D-Stars 

a ■ 
_ _ i 

To be specific, we shall work with a fixed model, where a global 0(3) x U em (l) symmetry is broken down to 
0(3) x U em (l). The Lagrangian density is (we work in units such that h = c = 1)[27], 

L = g^D^DvF + ^gTd^dvO* - if a a a a <jxj>* - ry 2 (00*) 2 - y (a a a a - a 2 ) 2 - ~ F^F^ (1) 

where = d^ + ieA^- F^ v is the electric-magnetic field vector, and F M „ = d^A^ — d^A^; a a is a triplet of scalar fields, 
isovector index a = 1, 2, 3; <f> is a complex scalar field with the U em (l) symmetry (f> — > e~ l6 '<fi. The field configuration 
describing a global monopole is 



o 
o 



cr a = cr f(p) — ,withx a x a = p 2 
P 



(2) 
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so that we will actually have a monopole solution if / — > 1 at spatial infinity. There is the current conservation 
3;^ = 0, where 



3" = ig^iF 



d<j> d(j>* 
d^ 7 ~ dx»' 



This leads to the conservation of the electric charge 



Introduce 



(3) 



(4) 



(5) 



so that / and h are both dimcnsionlcss and real. As the lowest energy solution of the theory, there should not exist 
electric current in and around the charged D-stars, so there should not exist the magnetic field. Therefore, the nonzero 
component of is A , which can be chosen as 



M{p) = -[v-g{p)] 
e 

The general static metric with spherical symmetry can be written as 



(6) 



ds 2 = B(p)dt 2 - A(p)dp 2 - P 2 (d0 2 + sin 2 Bdip 2 ) 



(7) 



with the usual relation between the spherical coordinates p,6,<p and the "Cartesian" coordinates x a . Introducing a 
dimensionless r = aop, from the Lagrangian density (1) and the definitions for /, h and g, the scalar field equations 
can be obtained as 



h + 



B_ A_ 

W ~ 2A 

2B 2A 



f = 



Af 
Ah 



-I + V 2 h 2 + 1 -X 2 (f 2 -1) 



(8) 



the equation for electric-magnetic field is 



.9 + 



B 
2B 



A_ 

2A 



g =Ae 2 h 2 g 



(9) 



where the prime denotes differentiation with respect to r, and A(r) and B(r) are the metric fields. By expressing the 
energy-momentum tensor of the system, one can reduce Einstein equation: G M1/ = SirGT^ . 

The charged D-stars consists of three regions: an interior surface and exterior, which can be discussed respectively 
[27]. For considering the exterior region of the star, one formally solves the Einstein equations of the static spherically 
symmetric metric as follows [28] 



A = 1 — 



8ttG 



Tip 1 dp 



(10) 



B 



A{ P ) 



exp 



SttG / (I? - TP)A{p)pdp 



where the time coordinate has been rescaled so that B = A 1 asp^oo. In terms of the dimensionless variable r 
and another dimensionless quantity e, the Einstein equation can be formally integrated and the solutions read as 



A- 1 (r) = l-e 2 - 



2Ga M A {r) , Ga 2 Q 2 A {r) 



+ 



Airr 2 
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B(r) = l-e 2 



2Ga M B (r) GalQ B (r) 



+ 



Airr 2 



(11) 



where 



M A (r) = 47rcr exp[-A(r)] f drexp[A(r)] x 

Jo 

[f-l + r 2 



+ ^ 2 (l- 2 + ^ i )(^ + / 2 )} 



, 2Ga M A , g" 2 



Q\ = 16^ 2 rexp[- A, (r)] J r 2 dr cxp [A x (r)](l - e 2 - — )j^ 



(12) 



and 



in which 



Mfl(r) = M A (r) exp[A(r)] + ^-^{1 - cxp[A(r)]} 



Q^ = Qiexp[A(r) 



A(r) = | jT dz(// 2 + /' 2 )z 



e 2 /" r , ff' 



AlW = 2 /„ 



A(r) - e 2 / <fe(/i' 

«/ OO 



2 , f '2 + 2A4- 5 2 /l 2 )^ 

2B ; 



(13) 



If this convergence is fast enough in asymptotic spacetime, M J 4(r)and Ms(r)will also quickly converge to finite values. 
Therefore, one can find the asymptotic expansions: 

/W = 1 -^-^+0(r- 6 ) 



M A (r) = M A + + 0(r- 3 ) 



(14) 



M B (r) = M A (r)(l - + + 0(r- 7 ) 

Since the effective mass M^(r) approaches very quickly its asymptotic value M^, it is a good approximation to 
take it as 



M A {r) — M A + 



(15) 



Therefore, to investigate the motion of light and a test particle around a charged D-star, we can take the metric 
coefficients as 



p p 



(16) 
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where, for convenience, we have rescaled M = u^Ma and Q = cr 2 a g7r — . 
3. The behavior of null geodesic outside the charged D-stars : 

The orbit of light outside of the charged D-star can be obtained by solving the geodesic equation 

d 2 x p p dx a dxl 3 _ 
a ? dr dr 



(17) 



. a j b 

where t is the affine parameter. But by using the fact that g a b^r^r = for null geodesies and the constant of 
motion 



and 



L = r 2± 



dr 

one can, considering the motion confined on the 9 = ^ plane, easily obtain the equation for geodesies as [29]. 

,L 2 



(18) 



(19) 



(20) 



Introducing the effective potential 



v l n(M L \ L 2 {l-e 2 ) L 2 GM GL 2 Q 2 

Veff = 2 B W(^) = ^3- + ~^A~ 



(21) 



the geodesies of light becomes the same as that of a neutral test particle with unit mass moving in the effective 
potential [30]. Taking G = 1 for convenience, it is easy to find that the effective potential has a maximum 



4L 2 (1 - e 2 ) 3 [3M 2 - 4Q 2 (1 - e 2 ) + M^9M 2 - 16Q 2 (1 - e 2 ) 



3M 



9M 2 - 16Q 2 (1 - e 2 ) 



(22) 



at 



r = r c = 



When 7}E 2 = V max , we will have 



3M+ v /9M 2 -16Q 2 (l-e 2 ) 
2(1 -e 2 ) 



(23) 



bcrit — 





3M+y / 9M 2 -16Q 2 (l-e 2 ) 


2 


2 \ 


/2(l-e 2 ) 3 


3M 2 - 4Q 2 (1 - e 2 ) + M^9M 2 - 16Q 2 (1 - e 2 ) 



(24) 



where b cr u is the critical value of 6, which is the generalization of the apparent impact parameter and is defined as 
b = jjj . The capture cross section of the charged D-star will be 

4 



7T 



07: 



irb 



3M + \J 9M 2 - 16Q 2 (1 - e 2 ) 



crit 



1 - e 2 ) 3 3M 2 - 4Q 2 (1 - e 2 ) + MJ9M 2 - 16Q 2 (1 - e 2 ) 



(25) 
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It is not difficult to prove that the path of light has a turning point at the largest radius, Ro, where ^\ r= n = 0. 
Using Eqs.(18)-(20), one can obtain the equation which can give the relation between R and b 

b~ 2 R% - (I - e 2 )R 2 + 2MR - 2Q 2 = (26) 
In order to compute the deflect angle of the light, we rewrite the Eqs.(18)-(20) as 

^ = [r%- 2 - (1 - e 2 )r 2 + 2Mr - 2Q 2 ]- 1 ' 2 (27) 
dr 

The change of light when passing a charged D-star should be Aip = tp^ — p-oo. Considering the symmetry, we have 

r°° dr 

Aip = 2 ^ = (28) 

Jro [r 4 b- 2 - (1 - e 2 )r 2 + 2Mr - 2Q 2 } 1 / 2 

The deflection of light up to the first order of M (M is small in the unit G = 1) is given by [29]: 

AM 



5tp = Atp-imM 1 ^' 



dM 



(1 - e 2 ) 3 / 2 i? 

M=0 



(29) 



4. The motion of timelike particles around a charged D-star 

From Eqs.(17)-(19), the equation of motion for a timelike test particle can be expressed as 

(~f=E 2 -?B(r)-^B(r) (30) 

where ju = ^- and E = ^- are the rescaled mass and energy of the test particle. Introducing \ = f ■> substituting it 
into Eq.(30) and then differentiating the equation with respect to ip, one will obtain the following equation: 

= i-(l- £ 2 + 7)x + «X 2 +/3x 3 (31) 
where p, 7, a and are dimensionless parameters defined as: 

1 GMJi 2 



(32) 



p L 2 
= V 2 GQ 2 

7 L2 J 7 " 
a = 3GM^ 

Noting that Eq.(31) is a nonlinear differential equation, it can be integrated formally as 

l fi- ( fo= / , = (33) 

J xo ^// x X o 2[I - (1 - 6 2 + 7 )x + a X 2 + 0x 3 ]d X + ±1 

where Xo i s the initial value of ^\ v=ipo . However, it is impossible to obtain the exact expression by integrating the 
above equation. In the following, we will gain some qualitative property of the system with the aid of phase-plane 
analysis without solving the equation numerically. To do so, we introduce two parameters as x = \ an d y — ^ and 
the autonomous system corresponding to Eq.(31) will be 
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% = f(*,y) = y (34) 

^ = g(x,y) = --(l-e 2 + 1 )x + a X 2 +(3x 3 
dip p 

Now, we prove the existence of periodic solution. According to the well known Bendixson's criterion [30], the 
equation of motion will have periodic solution if the divergence of the functional vector of the autonomous system is 
vanishing, i.e., V • (/, g) = 0. It is obvious that the functional vector (/, g) corresponding to the Eqs.(34) satisfies the 
criterion and therefore indicates that Eq.(31) has a periodic solution. One can also find that the periodic solution 
exists when e = 7 = (3 = which is the case of an ordinary star in asymptotically flat spacetime. This shows that the 
presence of constant and deficit angle and the electric charge will not exclude the existence of periodic solution from 
the equation of motion. 

Next, we analyze the critical points on the phase plane. The critical point is (xq, 0), where xq satisfies 

i - (1 - e 2 + + ax% + (ixl = (35) 

To analyze the type of the critical point, we firstly linearize the Eq.(34) and then do the translation x = x — xq. Thus 
the linearized equations should be: 

% - » < 36 » 

^ = Sx 
dip 

where 

6 = - (1 - e 2 + 7) + 2ax + 3f3x 2 (37) 
Using Eqs.(35) and (36), Eq.(37) could be rewritten as 

5= — — + 2(1 - e 2 + 7) - ax (38) 
pxo 

The eigenvalues corresponding to the system of equations will be 



Ai,2 = ± VS (39) 
The types of the critical point could be classified according to the eigenvalues as following: 

I. when S > 0, we have Ai < < A2, which indicates that the critical point is an unstable saddle point. Considering 
the Eq.(38), this case will correspond to the condition that 

(1) . a < and A < 0. 

(2) . a > 0, A > and (1^!±iWa < ^ < (1-^+7)+^ 

(3) . a < 0, A > and x > (l^dgWg > whcrc A = (1 - e 2 + 7 ) 2 - ^. 

II. when 5 < 0, we have two pure imaginary eigenvalues Ai,2 — ±iy/\S\, which indicates that the critical point is 
stable center. Considering Eq.(38), this case will correspond to the condition that 

(1) . a > and A < 0. 

(2) . a > 0, A > and x > ^^^ orO < x < , 

(3) . a < 0, A > and < x < (i-^+tO-Va ^ 

III. when S = 0, we have A1.2 = 0, which, together with the form of the autonomous system Eqs.(36), indicates that 
the motion is uniformly on the lines y = Constants and all the the points on the lines y = Constants are balanced 
positions. From Fig.l to Fig. 2, we show the phase graph for different initial values and different parameters. 



5. Discussion 
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FIG. 1: The phase graph when a — 0.100, 7 = 0.100, e = 0.010, p = 11 and 6 = 0.001, which corresponds to a charged case. 
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FIG. 2: The phase graph when a = 0.100, 7 = 0.000, e = 0.000, p = 11 and 8 = 0.000, which corresponds to an uncharged 
case. 

We study the motion of light and a neutral test particle around a charged D-star in this paper. We show that the 
deflect angle of light is different from the case in asymptotically flat spacetime for the existence of a factor of (1 — e 2 ) 
, which stands for the feature of a D-star that the spacetime has a deficit angle. 

Using the phase-plane analysis, we investigate the qualitative property of the dynamical equation controlling the 
motion of a test particle around the charged D-star. We prove that the equation of motion has a periodic solution, 
the position of whose critical point and type on the phase plane can be affected by the deficit angle and the charge. 
The conditions under which the critical point is stable center and unstable saddle point have also been given . 
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